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Abstract. MWC 560 (= V694 Mon) is the only known Symbiotic Star 
system in which the jet axis is practically parallel to the line of sight. 
Therefore this system is predestinated to study the dynamical evolution 
and the propagation of stellar jets. Spectroscopic monitoring done by 
Schmid et al. (2001) showed that the outflow is seen as absorption fea- 
tures in the continuum of the accretion disk and the accreting white 
dwarf, the emission line spectrum of the accretion disk and the spectrum 
of the red giant. We present the first numerical simulations of the jet of 
this particular object using the NIRVANA code (Ziegler & Yorke 1997) in 
order to reproduce the velocity structures seen in the observational data. 
This code solves the equations of hydrodynamics and was modified to 
calculate radiative losses due to non-equilibrium cooling by line-emission 
(Thiele 2000). 



1. The Observations 

The observational data revealed highly variable absorption structures which 
could be used to determine the velocity and acceleration structure of the jet. 
There exists a stable component with a radial velocity of —1100 =L 200 km s _1 
which represents the jet outflow at a larger distance. In addition, strong high 
velocity components with radial velocities up to —2500 km s _1 appear repeat- 
edly on time scales of a few days and decay within several days. At the same 
time, the absorption between the absorption trough and the emission line with 
radial velocities of ~ —900 to —400 km s" 1 disappears which is an indication for 
an acceleration process near the jet source (see dynamical spectra in Schmid et 
al. 2001). 



2. The Models 

Due to the enormous computational costs, only 2D-models with the assump- 
tion of axisymmetry (2.5D) have been calculated yet. Therefore the companion 
became physically a Red Giant "Ring". First, we started a purely hydrodynam- 
ical model without cooling. The parameters of this pulsed jet model are listed 
in Table 1. The influence of the stellar wind to the surrounding medium was 



1 



2 



Matthias Stute, Max Camenzind & Hans Martin Schmid 



simulated with an 1/r 2 law. The gravitational potential of the Red Giant and 
the Jet source was taken into account. The jet pulses were simulated by an 
enhancement of the jet velocity in the nozzle to 2000 km s _1 occurring every 
seventh day and lasting one day. Thereby also the jet density was varied to 0.25 
or 2 times the normal density. 



Table 1. Parameters of the models 



parameter 



value 

5 • 10 b cm" 3 
10000 K 
10 3 km s _1 
1.25 • 10 6 cm" 3 
1.0 • 10 7 cm" 3 
2 • 10 3 km s" 1 
1.0 • 10~ 6 M Q yr" 1 

10 km s -1 
1 AU 
1 AU 
4 AU 

50 x 30 AU 

20 grid cells/AU 



jet 



jet pulse 



stellar wind 



system parameters 



simulation parameters 



number density 

temperature 

velocity 

number density (models I & II) 
number density (model III) 
velocity (models I - III) 

red giant 
velocity 
jet-radius 

radius of the Red Giant 
binary separation 
integration domain 
numerical resolution 



In a second simulation, a cooling mechanism was implemented. According 
to Sutherland & Dopita (1993), three species (electrons, protons and hydrogen 
atoms) were considered and the heavier elements were simulated with a general 
cooling curve. In a third - again purely hydrodynamical - simulation, the num- 
ber density of the jet pulse was increased to eight times the value of the first 
two models. 



3. Results 

While the complete hydrodynamical simulation (450 days) took two weeks on a 
workstation, after nine weeks on a NEC SX5 vector computer only 74 days of 
the simulation with cooling were calculated. These computational expenses are 
caused by the unique parameter choice, especially the high densities, resulting 
in a small cooling time step. These are the first simulations in this parameter 
region including cooling. 

Comparing the logarithm of the total density for the same time step, one 
realizes the effects of cooling: due to smaller pressures in the jet the radial 
extension of the jet, the cross section and the resistance exerted by the external 
medium are lessened and therefore the cooled jet propagates faster (Fig. 1). 

Using the density and velocity data, theoretical line profiles can be calcu- 
lated, assuming the absorption I = Io e~ T for each grid cell with 
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T X = —\N j f jk . (1) 

m e c 



Figure 1. Logarithm of density after 74 days; without (left) and with 
cooling (right) 

With the increased number density of the jet pulse, the third simulation is 
able to reproduce the measured absorption of the fast pulse components (Fig. 
2). 

4. Discussion 

In the purely hydrodynamical model the high velocity component decays much 
faster than in the observations - on time scales of one day Taking into account 
the effects of cooling, it should decay slower which would be in agreement with 
the observational data. Due to the increasing propagation velocity of the jet and 
the decreasing resistance, the faster component should be able to exist longer. 
As this simulation is just at its beginning - and the jet has not yet reached 
its quasi-stationary phase, the results are still preliminary. But one first result 
could be the importance of cooling effects for modeling the jet correctly. 

A parameter study varying further the density and the injection velocity 
during the outbursts should bring deeper insights into the real mass loss rate 
and jet energy in this system. 
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Figure 2. Jet absorption line profiles for a series of eight consecu- 
tive days covering a high velocity pulse; theoretical (left) and observed 
(right) 



